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Abstract 

Inspired by the forms of delay-Painleve equations, we consider 
some new differential-discrete systems of KdV, inKdV and Sine-Gordon 
- type related by simple one way Miura transformations to classical 
ones. Using Hirota bilinear formalism we construct their new inte¬ 
grable discretizations, some of them having higher order. In partic¬ 
ular, by this procedure, we show that the integrable discretization of 
intermediate sine-Gordon equation is exactly lattice mKdV and also 
we find a bilinear form of the recently proposed lattice Tzitzeica equa¬ 
tion. Also the travelling wave reduction of these new lattice equations 
is studied and it is shown that all of them, including the higher or¬ 
der ones, can be integrated to Quispel-Roberts-Thomson (QRT) map¬ 
pings. 

pacs: 02.30.Ik, 02.30.Ks, 05.45.-a 


1 Introduction 

Integrable dynamics is a fascinating topic which has progressed and devel¬ 
oped rapidly in the last decades. Ronghly speaking, it is associated with 
the existence of a hnge amonnt of internal or hidden symmetry and this fact 
explains qnite well the predictability and regnlarity which characterize the 
integrable systems. Lattice eqnations (partial difference or PAE’s), which 
exhibit a higher complexity, attracted many stndies in the last years II1-0- 
The main progress was possible due to appearance of some efficient tools 
in proving integrability such as singularity conhnement ISE], ” cube consis¬ 
tency” [8]- [To] and complexity growth mm- 

Intensively studied in the past, the lattice versions of classical soliton 
equations have proven to be instrumental in establishing a classihcation of 
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quadrilateral lattice equations based on some symmetry assumptions (like 
cube consistency and tetrahedron property) [S] . However, in the case of higher 
order equations the situation is different. Because of the higher complexity 
they are studied mainly by examples, although in principle, from the Sato 
theory and Miwa transformations |15]-[T9], one can obtain lattice equations 
corresponding to a given hierachy. However there is another way to construct 
lattice equations starting directly from some differential-difference integrable 
versions and discretize directly by the Hirota bilinear method. It is easy 
to discretize the bilinear form but problems may appear at re-obtaining the 
nonlinear form. The integrability is proved by showing the existence of three 
soliton solution. The construction of 1-soliton and 2-sohton solution is easy 
in bilinear form and can also be done for some non-integrable systems. But 
the 3-soliton solution imposes strong restrictions on the equations and, ac¬ 
cordingly is practically equivalent with integrability. It was used succesfully 
for classifying integrable partial differential equations [25]-[28] and quite re¬ 
cently for PAE’s on a 3 X 3 stencil [32]. In this paper we are constructing 
some new lattice versions of KdV and mKdV equations starting from new 
differential-difference forms related to delay-Painleve equations. We use Hi¬ 
rota bilinear formalism for both construction and proving integrability. Also, 
starting from the bilinear form of semidiscrete sine-Gordon equation we find 
again the recently proposed lattice Tzitzeica equation. In the last section 
we make the travelling wave reduction of the above lattice equations and we 
show that all of them can be integrated to classical QRT mappings [31| 


2 New integrable lattice KdV equations 

We start from a new differential-difference equation namely: 

Wn = - Wn-l). (1) 

Wn—1 

The travelling wave reduction of this equation is exactly the autonomous 
limit of three point delay-Painleve I equation [20], [21]. There is also a one 
way Miura transformation from ([T]) to the well known differential-difference 
KdV equation of Hirota [21]; if Wn is a solution of ([I]), then = Wn/wn-i is 
a solution of 

hji rtfi—i) 0. (2) 
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Indeed if we put Un = Wn/wn-i, we obtain (after rearranging terms): 

^n—2{^n^n—l '^n^n+l) ^n(^n—l^n—2 l) 0- 

This relation is true since both brackets can be computed from ([1]) and 
its downshift along n. 

Integrability of ([T]) can be proved by Hirota bilinear formalism. If Wn = 
Fn+i{t)/Fn{t), then the following bilinear system is found (for notational 
simplicity we remove t-dependence): 


DtFn+i ■ Fn — Fn+2Fn-l — FnFn+l, (3) 

which is exactly the bilinear form of the classical integrable differential- 
difference KdV equation ([2]) [23]. This fact shows that the equation o 
is a potential semidiscrete KdV equation. 

In order to discretize a bilinear system we have to discretize the Hirota 
operator in (|2]) by replacing time derivatives with finite differences with the 
step 5 (and t 6m): 

^ F{n, 6m) — 


DtGn ■ Fn Q[G(n, (m 1)5) - G{n, 6m)] 


— ^^[F(n, (m -I- 1)5) — F{n, 6m)]^ G{n, 6m) 

After that, we impose the invariance of the resulting bilinear equations 
with respect to multiplication with exp(/in -|- um) for any p, u (the bilinear 
gauge invariance). 

The discretization of ([3]) is straitforward. We obtain the classical reduced 
Hirota-Miwa equation: 

rpm+l rpm _ rpm ipm+l _ i:t rpm+l rpm _ jpm rpm+l] 

^n+1 ^n+l-^n ~ n-1 ^n+1-^n 1' 


Consider = F)^i/F)^. Then, dividing the bilinear equation with 
we obtain the following: 


lym+l 


(l_5)Vp- + 5 




(4) 


which is a quad-lattice equation and it is the discrete analogue of ([I]). One 
can easily see that it is not consistent around the cube so it is not in the 
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ABS-list. To our knowledge, this is a new equation and, parallel to the 
semidiscrete case, there is a simple Miura transformation from it to classical 
lattice KdV equation of Hirota [23]. More precisely, if is a solution of 
dlj), then: 

< = (5) 

obeys the lattice KdV of Hirota: 


m+l _ym _ 


_ m+l _ 

“n \“n+l “ 


n—lJ 


( 6 ) 


Indeed, if -u™ = then (I6|) becomes: 


lO" [(1 - , + W”+>HC*+‘] - 


-Wn-2 [(1 - ] = o, 


which is true since the hrst square bracket can be computed from one down¬ 
shift along the n-direction of (jl]) and the second square bracket from the one 
upshift along the m-direction. 


3 New integrable lattice form of mKdV equa¬ 
tion 


In this section we are starting from the following new equation: 


Vn 


2,Vn 


"^n—l 


(7) 


which goes again in the travelling wave reduction to the three-point au¬ 
tonomous delay-Painleve II [20]. Also, there is a Miura from ([7]) to the 
classical semidiscrete mKdV equation (self-dual nonlinear network) [22112^ : 


hn (1 T ^n){,^'n+l '^n—l)i 


( 8 ) 


namely, = —logn^. This fact can be seen immediately since ([7]) is 
equivalent with: 


i d 

logUn = tan 
2dt ^ 


% % 

- log Vn+l - - log Vn-1 
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Applying tan“^ in both members and then taking the derivative with 
respect to t, one gets exactly ([H])- So the equation ([7]) is a potential new 
semidiscrete mKdV equation. 

We build the Hirota bilinear form using the substitution Vn = Gn{t)/Fn{t). 
Introducing in ([7]) and decoupling in the bilinear dispersion relation and the 
soliton-phase constraint, we obtain the following system: 

DtGn ■ Fn = Gn+lFn-l ~ G^-lF^+i, (9) 

‘^GnFn = Gn+lFn-l + Gn-lFn+l- ( 10 ) 


The above system is an integrable one since it admits 3-sohton solution 
of the following form (fcj is the wave number, uji is the angular frequency): 


Ail,A‘2,Ai3G{0,l} \i=l / i<j 

E (12) 

/il,/i2,Ai3G{0,l} \j=l / i<j 

where rji = kin + cuit, i = 1,3. 


The dispersion relation has the form: 

oji = (13) 

Phase factors and the interaction terms are: 


A. 


V 


hi — Q/i — 1, 

cosh(A:j — kj) — 1 
cosh(/i:j + kj) — 1 ’ 


f < j = 1,3. 


( 14 ) 

( 15 ) 


We discretize the bilinear system (|9l)- flT0ll in the same way. Replacing 
time derivatives in (Ej) with finite differences and imposing the bilinear gauge 
invariance, we obtain: 


G m+l 

n 


pm _ p,m p' 

J- ^ ‘Tt 'r 


m pm+1 
n n 




"Wn+l-^n-1 ^n-l-^n+1 J) 

^m rpm , /^m ip'm- 

'^n+l-^n-1 “T '-^n-l-^n+1- 


( 16 ) 

( 17 ) 
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The above system admits the following 3-sohton solution: 

Z (n“!“(pr9r")'“)n-4t'‘’. (18) 

A11,M2,A13S{0,1} \*=1 / i<j 

E (fl'>r‘(p??r‘)'“)npir' (19) 

Ml .M2,M3€{0,1} \i=l / i<j 


with the same phase factors and interaction terms as in the differential- 
difference case, but different dispersion relation: 

> = PO) 

where Pi = i = 1,2,3 (/c* is the wave number and cuj is the 

angular frequency). 

Now we can recover the nonlinear form. Dividing fITB]) by and 

using the following notations T™ = ”|i , we obtain: 

.m+1 m _ ^-pruf ,m+l m \ 


From equation (1T7|) we hnd: 

[^n ) ^ ^n-1 + ‘^n+1 

pm pm ~ 9/,,m 

But one can see immediately: 

^ FTA^F’n*' (Fnf _ 

Finally the nonlinear form of our sistem is: 


u 


m+1 

n 


-UJ. 


-pm 

n+1 


- K-i). 

< + <_! a.Si‘ j, 

<+2 + w,?+‘ wr 


( 21 ) 


One can eliminate F™ and F™_,_^ and we get the following new higher order 
nonlinear lattice mKdV equation: 


UJ. 


m+1 


UJJ. 


UJ. 


m+1 

n+1 


UJ. 


m+1 I m+1 

a)»,L0 UJn 


n-1 ^n+2 


UJ'. 


UJ. 


m+1 

n+1 


UJ. 


n+1 


UJ. 


m+1 

n+2 


UJ' 


/ } I ‘ 

OJ 


J- / 

+ ‘^n+l 


UJ. 


m+1 ■ 
n+1 


( 22 ) 
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4 Lattice system related to intermediate sine- 
Gordon equation 

In this section we are going to study the following differential-difference equa¬ 
tion: 

— {UnUn+l) = 'jul + (23) 

where 7 and k are constants. In the travelling wave reduction it goes to the 
bi-Riccati form of the delay-Painleve II equation [21]. On the other hand, for 
7 = —K = 1 the above equation is equivalent with the famous intermediate 
sine-Gordon equation ED: 

dtTy(x,t) + 2smy{x,t) = 0, 

where T is a singular integral operator defined through: 

(^/)(a^) = ^^y coth ^^^^^ f{z)dz. 

One of the key properties of the T operator is that it acts as a ’’jump” 
among different strips on the complex x-plane. So starting from a mero- 
morphic function u{x,t) in a horizontal strip of width 2a of the complex 
x-plane, we define y{x, t) = i \ogu/u where it = u{x -|- ia) and u = u(x — ia). 
Then Ty = —log mm, so the intermediate sine-Gordon becomes fIMl) for 
n = {x — ia)/2ia. 

The bilinear form of fl23ll can be computed by taking m„ = G„/ Fn. Taking 
into account the following relation (which can be checked by direct compu¬ 
tation). 


Gn+lFn+lDtGn’Fn+GnFnDtGn+l’Fn+l — Gn+lFnDtGn’Fn+l+GnFn+lDtGn+l'Fn 
the equation fl 2 ^ will turn into: 

Gn+lFn{DtGn ’ Rn+1 ~ f^Gn+lFn) + GnFn+l{DtGn+l ' Fn — yG^iG+i) = 0. 


Now it can be splitted in the following way: 


DfGn ■ Fn+l — nGn+lFn — AGnFn+i, (24) 

DtGn+l ■ Fn — 'yGnFn+l = —AGn+iFn, (25) 
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where A is a gauge constant. Replacing time derivatives with hnite differences 
and t with 5m in fl24p and 0251) . and imposing the bilinear gauge invariance, 
we contruct the fully discrete gauge invariant bilinear equations: 


or 1 - (26) 

U” - G”+iFr*‘ = '5(7G”JC+-i‘- 21G”«F„”)- (27) 


We take A=l, 7 = —k = 1 and the above system admits the following 
3-soliton solution: 


= E (n <■ ) n ^1 “'. 

/71,M2,/73S{0,1} \i=l / i<j 

Gr= E (n'>r(p”9rt'“)n7ir. 

/71,A‘2,/73G{0,1} Vi=l / i<j 


(28) 

(29) 


where the dispersion relation is: 


1 25 + Pi 

1 + (1 + 25)pi 


The phase factors are: 


ai = -bi = l, i = 1,3 


and the interaction terms: 


A,; 


f Pi-Pj y 

\i-PiPj) 


In order to see the nonlinear form, we take and from 

bilinear equation we obtain (after eliminating the term 


\rm _ 

^72+1 ~ 


j^m +1 


vrH-?(i+<^)-'^vr 

jf™(i + i) - 


which is nothing but classical lattice niKdV equation [29]. So, we can say 
that the lattice version of the polynomial intermediate sine-Gordon equation 
is the classical lattice mKdV. It is well known that lattice mKdV and lattice 
sine-Gordon are practically equivalent [33] (up to a exponentiation of the 





dependent variable with (—1)"*), so now we can say that lattice intermediate 
sine-Gordon is equivalent with lattice sine-Gordon (up to the above men¬ 
tioned exponentiation). 

Remark 1: 

Remark on discrete Tzitzeica equation: An interesting, although formal con¬ 
nection can be made with the discrete Tzitzeica equation which has been 
introduced quite recently by Adler [SU] . Starting from the general differential- 
difference equation (which goes to various forms of bi-Riccati delay-Painleve 
III in the travelling wave reduction): 

v! 11 

—-= auu H —- + Pou + fdiu (30) 

u u uu 

we get by means of u = G/F: 


DtG ■ G - fiFF = AGG, ( 31 ) 

DtF ■ F + aGG + ^qGF + ^^GF = AFF, ( 32 ) 

where A is a constant. 

We discretize in a gauge invariant way the above bilinear form and we 
get the following general bilinear system: 

^m+l^m /^m /^m+l r j^m+l T^m _ r A/^m /^m+1 /'QQ^ 

^n+1 ^n+l^n ^n+1 ~ 

TAm+1 T^m rpm ipm+l , ^m+1 , ro ^m+1 rpm , 

^n+l ^n+ln + + OPOO^n -^n+l + 

+i/3.„G™ 1^7+' + i/3„iG7F™y + ^/3„G7+‘F7 = (34) 

where (3oo, Poi, lo, h are arbitrary coefficients, which have to be determined 
according to integrability requirements. 

Now consider SA = = l,Sa = —l/c^,/Soo = —Poi = f/dc,(3io = 

—/dll = l/6c (this choice of coefficients will make more difficult the continuum 
limit (5 —> 0). Then the above bilinear system will have the form: 

^72+1 -^R+1 ~ 

„ Tpm+l jpm ^ G-im G-im+l i G-im+l ipm i rpm+l /^m rpm+l G-im+l rpm _ ri 

C-^n+l -^n ^n+l^n +^n -^n+1 + ^n+l-^n ^n-^n+1 ^n+l-^n ~ U. 

( 36 ) 
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Calling and eliminating the term we 

will get exactly the form of lattice Tzitzeica equation found by Adler |5U] : 

_ 1 _ 

c - - ^r^+y c-w™ - 'rr+‘ - irr+i 

We claim that the bilinear system fl35D . fl36l) is the Hirota bilinear form 
of the lattice Tzitzeica equation. Indeed, if we put G™ = = 

then the hrst bilinear equaton is identically verihed and the second 
bilinear equation is (up to a common factor) exactly the trilinear form found 
by Adler j30] : 


det 


m+2 _m+2 ^m+2 

' n '72+1 '72+2 

T-m+1 -1 m+1 T-m+1 

'n ^ 'n+l 'n+2 

-m 

'n 'n+l 'n+2 


(c 


,-l\ m m+l m+2 _ 
) n ' n+l ' n+2 


0 . 


5 Reductions of lattice equations 


5.1 Reduction of the potential lattice KdV 

In this section we are going to study various (p, g)-reductions or travelling 
wave reductions of the above lattice equations. We will show that even 
though some lattice equations are of higher order, their travelling wave re¬ 
ductions can be integrated to classical QRT-mappings [3l]. The interesting 
fact is that both additive and multiplicative QRT-mappings are obtained 
(additive and multiplicative refers to the type of rational surface on which 
the mappings act as automorphisms [3S])- The general form of a symmetric 
QRT mapping is the following: 


f2{Xm) - Xm-lfsiXm) 

where fi, f 2 , fs are general quartic polynomials in Xm- Any QRT mapping 
possesses an invariant which is biquadratic in Xm and Xm-i- The integrability 
comes from the fact that this biquadratic correspondence can be integrated 
in terms of elliptic functions (result which goes as far back as Euler). 

Let us start with the quad-lattice KdV equation ()4j): 


= (1 - 6)WT + "tR" . 

n—1 
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We consider that = x{n + m) = with u = n + m. In this rednction 
our equation becomes the following mapping: 

Xy+lXy-l - (1 - 5)XyXy-l “ 5Xy+2Xu+l = 0 . 


Dividing by x^j^iXy we get: 

Xy-l 


x„ 




X 1/+1 


X,, 


Wn—\ I Xu 


= l/Wy-l - (1 - 5)/{WyWy-l) 


5Wv+lWv 



Q = —1 + 1/5 


a + 1 a 

Wi,+iWy_l = - 


The last relation is a QRT mapping which is the autonomous limit of a 
g-Painleve equation realizing an automorphism of a rational surface of type 
R/)[35]. 

Remark2: The travelling wave reduction of classical lattice KdV of Hirota 
gives a QRT mapping which is the autonomous limit the d-Painleve equation 
(aditive type of rational surface -E'g^^)[35]. It was shown in [36] that, if = 
u{n + m) = Uy, 5' = 5/{I — 5), then we get from ([6|): 


1 1 


Uy+l + Uy + Uy-l — 7 ~ 77 - 

o' Uy 


What is really interesting is that the two equations (jl]) and ([6|) give, by 
the same reduction, both multiplicative and additive QRT mappings. 


5.2 Reduction of the higher order lattice mKdV 

The higher lattice niKdV is fl2^ : 


u 


m+l 

n 




m+l m 

n+1 ^n+1 


U! 


m+l 

n+1 


- CU, 


n—1 


U) 


m+l 

n+2 


+ 00 '. 


m+l 


U 


m 

n 


OO 


m+l 

n+2 


-UJl. 


OJ 


m 

n—1 




n+l 


OJ 


m+l * 
n+1 


We consider that u{n,m) = x{n + m) = Xy with u = n + m. In this 
reduction our equation becomes the following four-order mapping: 


Xy+l Xy Xy-\-2 Xy—\ Xy-^^ “j” Xy-^\ Xy 

Xy+2 Xy-^-l Xy-^-^ Xy Xy—l “t" Xy^l Xy-^-2 
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Then we have: 

-Xy _ Xy+2 - Xy-l Xy+^ + Xi,+ l X„ 

Xu-\-2 Xi/^i Xi,^^ Xi/ Xjy—i T Xi/_|_i Xj^_|_2 


_ Xu+l/Xu 1 Xu Xu^2lXu—\ 1 ^iz+s/^f+l T ^ Xu-t-1 

Xu+2/Xu+1 1 Xu-\-\ Xu-\-‘i/Xu 1 Xu-\-\jXu—I T 1 Xu-\-2 

_ Xu+i/Xu-1 Xu _ {Xu+2/Xu+l){Xu+i/Xu){Xu/Xu-i) - 1 
Xu+ 2 /Xu+1 - 1 Xu+1 {Xu+3/Xu+2){,Xu+2/Xu+l){Xu+l/Xu) - 1 

{Xu+3/Xu+2){Xu+2/Xu+i) + 1 Xu+i _ 

{Xu+l/Xu){Xu/Xu-l) + l Xu+2 

i»i/= 

Xi/ 

Wu - I _ Wu+lWuWu-1 - 1 Wu+2Wu+l + I Wu _ ^ 

Wu+1 - 1 Wu+2Wu+lWu - 1 WuWu-1 + 1 Wu+1 
Now we force {wuWu+i + 1 ) as a denominator in the second term. We get 


Wu-l _ Wu+lWuWu-l - 1 Wu+2Wu+l + 1 Wu+lWu + I Wu 

Wu+1 - 1 Wu+2Wu+lWu - 1 Wu+lWu + 1 WuWu-1 + 1 Wu+1 


, . {Wu - l){Wu+iWu + l){WuWu-i + 1 ) _ {Wu+i - l){Wu+2Wu+i + l){Wu+iWu + 1 ) 

{Wu+lWuWu-1 - l)Wu {Wu+2Wu+lWu - l)Wu+l 

Again the left hand side is the downshift of the right hand side, so both 
members are eqnal to a constant a. Accordingly, the integrated rednction of 
the lattice mKdV is the following more complicated QRT-form: 

I -Wu{(T+ 1) +Wu-l{Wu -wl) 

Remarks: Travelling wave rednction of Tzitzeica eqnation gives trivially 
the QRT mapping: 


Wu+lWu-l 


C Wi/^\ Wy — \ 

c~^wl — 2wu 
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6 Conclusions 


In this paper we have considered some differential-difference nonlinear equa¬ 
tions related to corresponding integrable nonlinear dynamical systems well 
known in mathematical physics like KdV, mKdV and sine-Gordon equations. 
For each of them, we have derived the Hirota bilinear form and introduced a 
Miura transformation that changes the new semidiscrete equations in already 
known differential-difference equations. Using Hirota bilinear formalism we 
have constructed their new integrable discretizations, some of them having 
higher order. The Hirota bilinear formalism was also used for proving integra- 
bility. The paper also shows that the integrable discretization of intermediate 
sine-Gordon equation is the lattice mKdV, and proposes a bilinear form of 
lattice Tzitzeica equation [30]. The travelling wave reduction of all these 
lattice equations leads to classical QRT mappings. 
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